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In this note, we will review and extend some results from our previous work [8] where we 
introduced a novel approach to imperfectly transferable utility and unobserved heterogeneity in 
tastes, based on a nonlinear generalization of the Bernstein-Schrodinger equation. We consider an 
assignment problem where agents from two distinct populations may form pairs, which generates 
utility to each agent. Utility may be transfered across partners, possibly with frictions. This general 
framework hence encompasses both the classic Non-Tranferable Utility model (NTU) of Gale and 
Shapley [6], sometimes called the “stable marriage problem”, where there exists no technology to 
allow transfers between matched partners; and the Transferable Utility (TU) model of Becker [T] 
and Shapley-Shubik m, a.k.a. “optimal assignment problem,” where utility (money) is additively 
transferable across partners. 

If the NTU assumption seems natural for many markets (including school choices), TU models 
are more appropriate in most settings where there can be bargaining (labour and marriage markets 
for example). However, even in those markets, there can be transfer frictions. For example, in 
marriage markets, the transfers between partners might take the form of favor exchange (rather 
than cash), and the cost of a favor to one partner may not exactly equal the benefit to the other. 

In [8], we thus developped a general Imperfectly Transferable Utility model with unobserved 
heterogeneity, which includes as special cases the classic fully- and non-transferable utility models, 
but also extends to collective models, and settings with taxes on transfers, deadweight losses, and 
risk aversion. As we argue in the present note, the models we consider in [8] obey a particularly 
simple system of equations we dubbed “Nonlinear Bernstein-Schrodinger equation”. The present 
contribution present a general result for the latter equation, and also derives several consequences. 
The main result is derived in Section I. Section 2 and Section 3 consider equilibrium assignment 
problems with and without heterogeneity. Finally, we provide a discussion of our results in Section 
4. 

* Economics Department, Scienees Po, Paris, and New York University. Email: alfred.galichon@sciencespo.fr. 
Galichon gratefully aeknowledges funding from the European Research Council under the European Union’s Seventh 
Eramework Programme (FP7/2007-2013) / ERC grant agreements no 313699 and 295298, and FiME. 

^Society of Fellows, Harvard University. Email: kominers@fas.harvard.edu. Kominers gratefully acknowledges the 
support of NSF grant CCF-1216095 and the Harvard Milton Fund. 

^Economies Department, Scienees Po, Paris. Address: Seiences Po, Department of Economies, 27 rue Saint- 
Guillaume, 75007 Paris, Franee. Email: simon.weber@scieneespo.fr. 


1 



1 The main result 


For X G and y € y, we consider a function Mxy : ^ M, and Mxo : M ^ M and M^y : M —>■ 

M. Let {nx)x;fzx {'i^y)y^y be vectors of positive numbers. Consider the nonlinear Bernstein- 
Schrodinger system, which consists in looking for two vectors u € and v € such that 

f y^xO i^x) “ 1 “ ^^y^y y^xy {Ux, Vy) — Tlx 
1 yiOy {Vy) + 'Yhx&X ^xy {Ux,Vy) = ITly 

Theorem 1. Assume Mxy satisfies the following three eonditions: 

(i) Continuity. The maps Mxy : {ux,Vy) i— Mxy{ux,Vy), Mxo : {ux) '— Mxo{ux) and Moy : 

(vy) I— yioy{vy) 0 X 6 conUnuous 

(a) Monotonicity. The map Mxy : {ux,Vy) i— Mxy{ux,Vy) is monotonically decreasing, i.e if 
Ux < u'x and Uy < u'y, then Mxy{ux,Vy) > M{u'x,,v'y). The maps Mxo : (ux) '— Mxo{ux) and 
Moy : (vy) I —> MQy{vy) are monotonically decreasing. 

(Hi) Limits. For each Vy,lim.u„,^oo Mxy{ux,Vy) = 0 and Mxy{ux, Vy) = +oo, and for 

eachux, liuiuy^oo Mxy{ux,Vy) = 0 and liuiuy^-oo Mxy {ux,Vy) = +oo. Additionally, Mxo{ux) = 

0 and liuiuy^oo Moy{uy) = 0 , and Mxo{ux) = +oo and liuiuy^-oo Moy{uy) = +oo. 

Then there exists a solution to the nonlinear Bernstein-Sehrodinger system. Further, if the 
maps M are , the solution is unique. 


This theorem appears in [8] under a slightly different form. The proof is interesting as it 
provides an algorithm of determination of u and v. It is an important generalization of the Iterated 
Projection Fitting Procedure (see and m), which has been rediscovered and utilized many 

times under different names for various applied purposes: “RAS algorithm” m, “biproportional 
fitting”, “Sinkhorn Scaling” [3], etc. However, all these techniques and their variants can be recast 
as particular cases of the method described in the proof of Theorem [TJ For convenience, we recall 
the algorithm used to provide a constructive proof of existence. 


Algorithm 1. 

Step 0 
Step 2t + 1 

Step 2t + 2 


Fix the initial value of Vy at Vy = +oo. 

Keep the values Vy^ fixed. For each x € T, solve for the value, such 

that equality J2yey ^xy{ux,Vy^) + Mxo{ux) = rix holds. 

Keep the values fixed. For each y € y, solve for which is the value, 

Vy^~^‘^ such that equality YIxgx ,Vy) + MQy{vy) = my holds. 


Then 


w and u* converge monotonically to a solution of the Bernstein-Sehrodinger system. 


In practice a precision level e > 0 will be chosen, and the algorithm will terminate when 
supj^ \Vy^~^‘^ - Vy^\ < e. 

Proof, (i) Existence. The proof of existence is an application of Tarski’s fixed point theorem and 
relies on the previous Algorithm. We need to prove that the construction of and v‘^~^‘^ at 

each step is well defined. Consider step 2t + 1. For each x G A, the equation to solve is 


Mxy{Ux, Vy) + MxoiUx) = Ux 

y&y 


but the right handside is a continuous and decreasing function of Ux, tends to 0 when Ux —>■ +oo 
and tends to +oo when Ux —)• —oo. Note that by letting Vy —>■ +oo, the terms in the sum tends 
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to 0, providing a lower bound for Ux- Hence is well defined and belongs in (^M^Q{nx),+oo^, 

and let us denote 

= Fxiv^^) 

and clearly, F is anti-isotone, meaning that for all y G implies Fx{v‘^^) < Fx{v^^) for 

all X € X. By the same token, at step 2t + 2, is well defined in (^M^y^^niy), -|-oo^, and let us 

denote 

where, similarly, G is anti-isotone. Thus 


„2t+2 = CoF (n2*) 

where G o F is isotone. But Vy<oo = Vy implies that < G o F Hence is a 

decreasing sequence, bounded from below by 0. As a result converges. Letting n, its limit, 

and letting u = F{v), one can see that {u,v) is a solution to the nonlinear Bernstein-Schrodinger 
system, (ii) Unicity. Introduce map C defined by 


f Cx — FIxy {Ux, Vy') F Mxo{ux)\ 

VCy = YlxeX F^xy {Ux, Vy) + FlQy{Vy)) 


One has 


where: 


F(^ {Ux, Vy) - 


A 

B 

G 

D 


• A = (dCx/dux')xx' = Yjy'&ydu^FIxy' {ux,Vy') +l\ix = x',t) otherwise. 

* B i^Cx / dVy) ^y dxyFIxy {UxjVy) 

* G id(^y/dUx)y^ du^FIxy {UxjVy) 

• D = {dC,y/dvy>)yy, = Y^x'&x dvyFIx'y {ux'Vy) + liiy = y',Q otherwise. 

It is straightforward to show that the matrix DC, is dominant diagonal. A result from m states 
that a dominant diagonal matrix with positive diagonal entries is a P-matrix. Hence DC {ux,Vy) is 
a P-matrix. Applying Theorem 4 in [9] it follows that C is injective. □ 


2 Equilibrium Assignment Problem 

In this section, we consider the equilibrium assignment problem, which is a far-reaching generaliza¬ 
tion of the optimal assignment problem. To describe this framework, consider two finite populations 
I and CT, and a two-sided matching framework (for simplicity, we will call “men” and “women” 
the two sides of this market) with imperfect transfers and without heterogeneity. Agents i € T and 
j G CT get respectively utility Ui and Vj they get at equiliburium. If i or j remains unmatched, they 
get utility 0; however, if they match together, they may get any respective utilities Ui and Vj such 
that the feasibility constraint is imposed 


'^ijiui,Vj) <0, ( 1 ) 

where the transfer function 'I'jj is assumed to be continuous and isotone with respect to its ar¬ 
guments. Note that at equilibrium, rij > 0 and Vj > 0 as the agents always have the option to 
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remain unassigned; by the same token, if for any pair i, j (matched or not), one cannot have a strict 
inequality in ([T]), otherwise i and j would have an incentive to form a blocking pair, and achieve a 
higher payoff than their equilibrium payoff. Thus the stability condition ^ij (ui,Vj) > 0, holds in 
general. Let fiij = 1 if i and j are matched, and 0 otherwise; we have therefore that Hij > 0 implies 
that Tjj {ui,Vj) = 0. This allows us to define an equilibrium outcome. 

Definition 1. The equilibrium assignment problem defined by T has an equilibrium outeome 
{^ij,Ui,Vj) whenever the following conditions are met: (i) Hij > 0, rt* > 0 andvj > 0 (ii) YljTij ^ 1 
and YliTij ^ '^ij ^ 0 (iv) fiij > 0 implies Tjj {ui,Vj) = 0. 

Note that, by the Birkhoff-von Neumann theorem, the existence of an equilibrium in this 
problem leads to the existence of an equilibrium satisfying the stronger integrality requirement 
Hij G {0,1}. This general framework allow us to express the optimal assignment problem (match¬ 


ing with Transferable Utility), as the case where 


Tjj {ui, Vj) = Ui + Vj - ^ij, 

while in the NTU case 

{ui, Vj) = max(Mi - aij,Vj - 'jij). 

Other interesting cases are considered in [8j. For instance, the Linearly Transferable Utility {LTU) 
model, where 

'kij (iti) Vj) — \ij{Ui ®ij) T Cijivj T*i) 


with Xij,Cij > 0) and the Exponentially Transferable Utility (ETU) model, in which Tjj takes the 
form 


'^ij (ui,Vj) = rlog 


exp(ni/r) -|- exp(uj/r) 


In the ETU model, the parameter Tij is defined as the degree of transferability, since r —>■ -t-oo 
recovers the TU case and r ^ 0 recovers the NTU framework. 


Theorem 2. Assume T is such that: (a) For any x £ At and y £ y, we have Afxy (•, •) continuous, 
(b) For any x£X,y£y,t<t' and r < r', we have 'I'j,y(t,r) < xy {f ,r'); furthermore, 
when t < t' and r < r', we have ^xyit^v) < A>xy{t',r'). (c) For any sequence {tn,rn), if (vn) is 
bounded and tn +oo, then lim inf Atxy (tn, Vn) > 0. Analogously, if (tn) is bounded and r^ -|-oo, 
then liminf (f„, Tn) >0. (d) For any sequence {tn,rn) such that if {tn — Vn) is bounded and 
tn —^ —oo (or equivalently, —>■ -|-ooj, we have that lim sup {tn,rn) < 0. Then the equilibrium 
assignment problem defined by Tjj has an equilibrium outcome. 

Proof. Consider T > 0 and let 


Mij {ui,Vj) = exp 
Mio {ui) = exp 
Moj (vj) = exp 

and consider the Bernstein-Schrodinger system 


Atij {Ui,Vj] 


-(fl) 

t) 


( Mio (ui) + Y(j&j ^ij '>Jj) = 1 
\ Moj (Vj) + 1 
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We need to show that Mxy {.,.), Mxo{.) and Moy(.) satisfy the properties stated in Theorem[TJ It is 
straightforward to show that conditions (i) and (ii) in Theorem [T] follow directly from assumptions 
(a) and (b) on T. Moreover, letting T ^ 0~^, assumptions (c) and (d) together imply that condition 
(iii) is satisfied. Hence we can apply Theorem [H and it follows that a solution uf, vj to the system 

exists. Note that Mjo {uf) < 1 and M^j < 1 imply that uj > 0 and vJ > 0. Now, consider 

the sequence obtained by taking T = k, k £ N. Then, up to a subsequence extraction, we may 
assume ^ Ui £ M’*' U {+oo} and Vj £ M"*" U {+oo}. It follows that converges 

in M"'' U{+oo}, hence = Mij converges toward fiij £ [0,1]. Similarly, the limits Jlio and 

floj exist in [0,1]. Hence (i) in Definition [T] is met. By continuity, ft satisfies 

f P-iO + — 1 

I Poj + Pij ~ ^ 

which established (ii). Let us show that (iii) holds, that is, that {ui,Vj) > 0 for any i and j. 
Assume otherwise. Then there exists e > 0 such that for k large enough, Tjj (u^, < —e, so that 

> exp(e/r) —>■ + 00 , contradicting fiij < 1. Thus, we have established (iii). Finally, we show 
that (iv) holds. Assume otherwise. Then there is i and j such that fiij > 0 and {ui,Vj) > 0. 

This implies that there exists e > 0 such that for k large enough, > e, thus Tjj (u^,Vj'^ < 
—Tloge —)■ 0, hence 'k^- {ui,Vj) < 0, a contradiction. Hence (iv) holds; this completes the proof 
and establishes that (/2, u, v) is an equilibrium assignment. □ 

3 ITU matching with heterogeneity 

Following [8], we now assume that individuals may be gathered into groups of agents of similar 
observable characteristics, or types, but heterogeneous tastes. We let A and y be the sets of types 
of men and women. An individual man i £ X has type Xj € A; similarly, an individual woman 
j £ J has type pj £ y. We assume that there is a mass Ux of men of type x and niy of women of 
type y, respectively. Assume further that 

'kij {Ui, Vj') — '^Xiyj {yi Vj Xfjxj') i 

where e and rj are i.i.d. random vectors drawn from a Gumbel distribution, and where T > 0 is 
a temperature parameter. Unassigned agents get Tejo and Tijoj. For all i such that Xi = x and 
Pj = y, the stability condition implies 

{ui X£iy, Vj Xrixj) ^ 0. 

Hence, 

min 'iixiVi {ui - Tsiy.Vj - Trjxj) > 0. 

i\Xi=x 

3-yj=y 

Thus, letting 

Uxy = min {ui -Tsiy} and Vxy = min {vj - Tpxj} , 
i:xi=x j-yj=y 

we have ^^xy (Uxy, Vxy) > 0, and with 

k'xy = l^ij 

i:Xi=X 

j-yj=y 
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we have that fi^y > 0 implies '^xy (Uxy, Vxy) = 0, and by a standard argument (the random vectors 
e and r] are drawn from distributions with full support, hence there will be at least a man i of type 
X and a woman j of type y such that i prefers type y and j prefers type x, that is, y-xy > 0 for all 
X and y) 

^ xy (Uxy, ^xy ) = 

Note that this is an extension to the ITU case of the analysis in Galichon and Salanie [7], building 
on Choo and Siow [3]. We have that 


Ui = max {Uxy + T£iy,T£io} and vj = max{I4y + Tr]xj,Trioj} , 


thus a standard result from Extreme Value Theory (see Choo and Siow [3] for a derivation) yields 


so we see that fixy satisfies 


Uxy = riog ^ and Vxy = riog 

y-xO IJ-Oy 


f^xO T f^xy — 


MOy + YIxGX — 


m„ 




xy 


(riogi^Tlog!^) = 0 
\ ^ ^ yoy J 


( 2 ) 


The various cases of interest discussed above, namely TU, NTU, LTU, and ETU cases yield, re¬ 
spectively, 

1/2 1/2 ^xy /rpTTX 

yxy = yJo yoy exp (TU) 
yxy = min {yxoe°‘^y , yoye ^^^) (NTU) 

^ _ ^VxyOLxy+C,xylxy)/Vxy+<^xy) ^^^y/V^y+Cxy) ^Qxy/Vxy+Cxy) 

( g,-axylTx:y _j_ Q-^xy j Tx,y V^y \ 


f^xy — 


(ETU), 


(see my To apply Theorem [U we let Mxy {ux,Vy) be the value m solution to (for a proof of 
existence and uniqueness of such a solution, see Lemma 1 of [8]) 


^xy {TlogniV Ux,Tlogm + Vy) = 0, 


and let 

Mxo (ux) = exp and Moy (vy) = exp 

In [8], we rewrote system (l2|) as a nonlinear Bernstein-Schrodinger system, namely 

^xO {Ux) T ^xy {Uxj '^y') ~ ^x 

Moy (Vy) + YIxGX ^xy {Ux,Vy) = TTly 

Theorem 3. The nonlinear Bernstein-Schrodinger system in ^ has a unique solution 

Proof. The proof follows directly from the application of Theorem [TJ It is easy to check that the 
conditions on Mxo{.) and Moy(.) required by Theorem [1] are met in this case. Lemma 1 in [8] 
provides a proof that Mxy satisfies these conditions. □ 
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4 Discussion 


In this note, we have argued how matching problems may be formulated as a system of nonlinear 
equations, also known as the Bernstein-Schrodinger equation or Schrodinger’s problem [15]. We 
have shown existence and uniqueness of a solution under certain conditions, and have explicited 
the link with various matching problems, with or without heterogeneity. Solving such a system 
of equations requires an algorithm that we call the Iterative Projection Fitting Procedure (IPFP); 
in practice, this algorithm converges very quickly. Our setting can be extended in several ways. 
One of them is to consider the case with unassigned agents. In that case, we have the additional 
constraint that ~ thus the nonlinear Bernstein-Schrodinger system, which in this 

case writes as 

f Xyj/ey Vy) = Ux 

1 '^xex ^xy {ux^Vy) = my 

has a degree of freedom, as the sum over x G of the first set of equations coincides with the sum 
over y (z y ol the second one. The one-dimensional manifold of solutions of this problem is studied 

in[2]. 
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